Single-Spin Addressing Using a Spatially- Varying Dressing Field 
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We propose a simple method to spectrally resolve single-spins in a cold atomic system, thus real- 
izing single-spin addressing. This scheme uses a dressing field with a spatially-dependent coupling 
to the atoms. We realize this scheme experimentally using a linear chain of trapped ions that are 
separated by ^ 3 /xm, dressed by a laser field that is resonant with the micromotion sideband of a 
narrow optical transition. 
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In recent years, remarkable capabilities of control and 
probing of atomic systems have been achieved. Ions and 
atoms are now routinely imaged at the single-particle 
level, enabling the study, among other topics, of many- 
body physics at an unprecedented level of accuracy [H,[2|. 
For many such applications, single-particle addressing 
is an increasingly important tool. Early realizations of 
individual-atom addressing relied on large inter-atomic 
separations 0, HI, However, applying single-addressing 
schemes for short inter-particle spacings is very challeng- 
ing, yet often a necessary task. For example, in order 
to investigate strongly interacting quantum systems, 
one has to reduce inter-particle spacings down to near- 
diffraction-limited distances. Significant experimental 
efforts have recently been invested in order to develop 
robust individual-particle addressing schemes with high 
spatial resolution, both in trapped neutral and ionic 
systems. Most setups rely on two general approaches: 
(i) A laser is tightly focused so as to interact only with a 
single particle, either directly manipulating its quantum 
state [4|, |5[), or inducing a light-shift which spectrally 
distinguishes the frequency of that atom from the rest 
(ii) In analogy with magnetic resonance imaging 
techniques, a strong gradient in magnetic field is applied 
in order to Zeeman-shift the atomic levels of interest 
depending on the position of the particle. Single atoms 
can be addressed using the transition frequency that 
matches their location [71-llOj. While the first approach 
is very demanding in term of laser spatial stability and 
resolution, and becomes difficult when inter-particle 
spacing compares with the optical wavelength of the 
addressing setup, the second requires the generation of 
sufficiently strong magnetic field gradients in order to 
resolve the frequency separation between transitions in 
closely-separated particles. 

In this letter we propose a simple scheme for single 
spin addressing in a cold atomic ensemble, that relies 
on a spatially-dependent dressing field. The resulting 
dressed eigenstates have spatially-dependent energies 
and therefore, transitions between dressed states in 
atoms at different locations are spectrally resolved. We 
realized this scheme using a linear array of trapped ions. 
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FIG. 1: (Color online) The system considered is an array of 
N identical non-interacting two-level systems, subjected to an 
external coupling field with a Rabi frequency Qj for particle j. 
The lower panel shows a sketch of the single-spin addressing 
protocol, involving three steps: (i) the bare eigenstate of all 
particles t and ^ at the Bloch sphere poles are mapped onto 
the dressed eigenstates in the sphere equatorial plane (shown 
by thick black arrows), (ii) A driving field (shown by the 
green arrow) which is amplitude-modulated at a frequency Qj 
rotates the spin of particle j (shown by the thin red arrow, 
and the red two- level system in the upper panel), (iii) The 
dressed states are mapped back on the bare basis at the sphere 
poles. All particles end up in state |, except particle j which 
is in t- 



The dressing field was provided by the inhomogenous 
micromotion of ions in the chain and a narrow linewidth 
laser that was tuned to the micromotion sideband 
of a narrow optical quadrupole transition. Quantum 
information manipulation using microwave-dressed ion 
qubits has been previously demonstrated in 

Different techniques for using the spatial variation of 
micromotion in ion traps for the purpose of individ- 
ual addressing have been suggested and demonstrated 
Here, we consider a non-interacting array of N 
two- level systems that are held at fixed positions. The 



Hamiltonian for particle j is H^j^ = 
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<j z the z-Pauli matrix, and H. 
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the interaction Hamiltonian between the atom and the 
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laser field, = ftVtj cos(kxj — ujLt)<jy , where uol 

is the frequency of the laser field, and /c, its wavenum- 
ber. Note that while the transition energy Tiujq is as- 
sumed equal for all particles, the Rabi frequency associ- 
ated with this coupling Qj is particle-dependent. In the 
rotating-wave approximation (RWA) and assuming the 
Lamb-Dicke regime, the Hamiltonian in the interaction 
picture reads (l4j . 

= y- ™h u) a) 

^dress 2 y ' [ } 
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where the frequency of the laser is tuned to the atomic 
transition uj l = ujq. In order to individually address par- 
ticle n, we note that is formally equivalent to the 
' dress J H 

Hamiltonian of an array of spin- 1/2 in a fictitious mag- 
netic field gradient, in the ?/-basis. In other words, the 
two-level systems are coupled to a spatially varying dress- 
ing field, and the energies of the dressed states, which are 
(i) 

the eigenstates of a y , are position-dependent. 

This analogy enables to construct a simple protocol 
for individual spin rotations, as is illustrated in FigHJ (i) 
The bare spin eigenstates (a z basis) are mapped onto 
an eigenstate of the dressed spin (<j y basis) |t)d/|i)d? 
and the strong resonant dressing field is turned on, 
locking the eigenstates in place by opening an energy 
gap between them. This energy gap is determined by 
the local coupling strength of each atom to the dress- 
ing field, (ii) A small <j x perturbation is amplitude- 
modulated (in the c^o-rotating frame), with a modu- 
lation frequency 5, so that the Hamiltonian becomes 

^dress = ^dress + EjLi ^- cos(<5t)^'>. This Hamil- 
tonian is solved for each particle n by performing an 
additional change to a rotating frame, at an angular fre- 
quency £l n , assuming Aj, \S — Qj\ <C <S, and neglecting the 
counter-rotating terms in this frame of reference [15|, [L6| . 
In that case, we readily find that the probability for spin 
j to be I t)dj while all the others are | \)d is 

p(t)=p J (t)H(i-p l (t)), (2) 

where Pi(t) = ((A^/4^) 2 sin 2 (^t/2) is the Rabi preces- 
sion, at a frequency = yAf -f A 2 /4 and the detuning 
in the ^-rotating frame is A$ = 5 — Qi. For all fij's 
different from (5, the spins remain locked on y-ax.is while 
the spin of particle n is resonantly rotated around the 
x-axis. (iii) The dressed states are mapped back to the 
a z bare basis. Following this procedure, all particles 
remain in their initial state except for particle n, whose 
spin has been rotated. 

We experimentally demonstrated the above individual 
spin-addressing scheme using a chain of atomic ions, 
held in a linear Paul trap, and typically separated by 
~ 3 jam. The experimental setup is presented in detail 



elsewhere [ItJ- In short, we trapped and Doppler-cooled 
one or more 88 Sr + ions in a linear Paul trap. We 
implemented a two- level system using the 5Si/2,+i/2 
state, hereafter referred to as the bright state, due 
its photon statistics during state-selective fluorescence 
detection [lij ]. and the long-lived 4D 5 / 2 ,+3/2 state, here- 
after called dark state. These two states are connected 
by an optical, electric-quadrupole, transition at 674 
nm. This transition is driven using a narrow-linewidth 
diode laser stabilized to an external high-finesse cavity 
[18|. In an ideal linear Paul trap, the micromotion, 
due to the ac-field modulated on the RF electrodes, 
vanishes along the symmetry axis of the trap. However, 
rather than a line of nulled RF amplitude as in the 
ideal case, the boundary conditions imposed by the trap 
end-caps generate a linearly increasing RF amplitude 
around a null- amplitude point at the trap center [l9| . 
The inhomogenity of this excess micromotion along the 
chain is used as a spatially-dependent coupling of the 
ions to the laser field. When operating the laser on 
the micromotion sideband of the quadrupole transition, 
located at a detuning of ^j^p = 21.75 MHz from the 
carrier transition, the resulting Rabi frequency depends 
on the ion micromotion amplitude via: Qj = Q c Ji(Vj)i 
where Q c is the carrier Rabi frequency, in the absence 
of modulation, J\ is the first-order Bessel function 
of the first kind [20]. The parameter r]j = k • Xj is 
the micromotion Lambe-Dicke parameter of particle j 
along the laser wavevector k, and Xj is its micromotion 
amplitude. For micromotion amplitudes small compared 
with the transition wavelength, rjj <C 1, the local Rabi 
frequency can be approximated by Qj ~ r)j.Q c /2. The 
micromotion amplitude of the ions in the chain can 
be controlled by applying a differential voltage to the 
trap end-caps and displacing the ions along the trap 
axis, or by varying the RF voltage on the trap electrodes. 

The values of the dressed state energies TiVtj for the 
different ions was first calibrated. The strength of the 
dressing field was measured through the frequency of the 
Rabi oscillations on the micromotion sideband for differ- 
ent end-cap voltages. For multi-spin systems, the popu- 
lation of ions in the dark state was estimated from photon 
detection histograms. FigEk shows the fraction of ions 
in the dark state versus the micromotion sideband pulse 
time, for a three-ion chain. The different Rabi frequen- 
cies involved can be conveniently extracted by Fourier 
analysis of the Rabi nutation curve, as shown in Fig|2]3. 
Three frequencies, differing by roughly 5 kHz are clearly 
observed. In Fig|2t, such spectra are shown as columns 
as a function of the global axial displacement of the chain. 
As seen, the three frequencies shift linearly with the dis- 
placement. This shows negligible effect of the spatial in- 
homogeneity of the laser beam profile on the dressed state 
energies on the length scales of Fig|2] One frequency nulls 
every time one ion sits on the RF-null (points labeled 1 



3 




30- 



(c) 



ffi 20 



g 10 
fa 



10 20 30 40 50 60 70 
Frequency (kHz) 



k0 



2 4 6 8 10 12 14 
Displacement (yum) 




1 2 3 4 5 

Addressing Pulse Phase (rad) 



FIG. 2: (Color online) Calibration of the single-spin address- 
ing protocol, (a) and (b) show the time evolution, and Fourier 
spectrum respectively of the fraction of dark ions in the three- 
ion chain driven on the micromotion sideband, for a given ax- 
ial position of the chain, (c) Calibration of the dressed state 
energies by Fourier transform of the Rabi nutation curve on 
the micromotion sideband as a function of the axial displace- 
ment. The three particular points on the map (labelled 1 to 
3) correspond to the RF-null being respectively on the central 
ion, in between two right-most ions and on the right-most ion, 
as shown in the schematic drawing on top of the map. The 
ions (RF-null) are represented by red (white) disks, (d) 2D 
scan of the phase of the amplitude-modulated driving field 
(with respect to the dressing beam phase) and the pulse time 
for a single ion. A a x modulation is obtained for a phase 
cj) ~ 3.03 rad (vertical dashed black line). 



and 3 on FigGfc). 

In order to apply this scheme to a multi-spin system, 
one has to address simultaneously all spins for the steps 
(i) and (iii). This can be solved in two different ways. 
(1) If micromotion is sufficiently weak so that the carrier 
Rabi frequency is almost constant throughout the chain, 
then one can apply a collective carrier 7r/2-pulse (along 
a x ) to bring the spins to the equatorial plane. Here, 
the carrier and the micromotion sideband transitions, 
together with the RF signal for the trap electrodes 
must all be phase-locked. (2) The spins are set on 
the equatorial plane using Rapid Adiabatic Passage 
(RAP) on the micromotion sideband interrupted on 
resonance. The RAP beam is kept on resonance, and 
thus automatically dresses the spins. In both schemes, 
following initialization, a AM-a^ pulse rotates the spin 
of the particle of interest. After the single-spin rotation 
the dressed states are mapped back to the bare states 



by reversing step (i). While we have successfully imple- 
mented both methods, we preferred the former, since 
RAP is several times slower than the carrier 7r/2-pulses. 

The amplitude-modulated cr x field was generated by 
adding a radio-frequency, that was amplitude-modulated 
at a frequency 5, to the laser acousto-optic modulator 
(AOM). This signal had the same carrier frequency 
and is phase-shifted by tt/2 from the dressing signal. 
Since the cr^-dressing and the cr x driving field were 
provided by the same beam, then = a^i, and the 
Rabi frequency is proportional to Qi when the 
detuning vanishes. The phase of the driving field 
was tuned in the following way: after the spins are 
dressed, the phase and time of the modulation are 
scanned, while the driving frequency 5 is tuned to a 
single-spin resonance. The resulting spin-flip probability 
for a single ion is displayed in Figj2]i. The period of 
the pattern is tt and we observed that at phases ft and 
ft + 7r (where ft ~ 1.46) the scan is independent of 
time, corresponding to a modulation colinear with the 
dressing field (</> = dzy). At phases ft ± ir/2 (vertical 
dashed line in Figj2]i), we observe the highest contrast 
for the oscillations, corresponding to resonant a x rota- 
tions of the dressed state. Note that the sum of the 
cr^-dressing and the cr^-modulation could be equally well 
produced by direct frequency modulation of the a y beam. 

We next scanned the frequency of the driving field 400 
/is pulse. A typical spectrum, displayed in Figj3^, ex- 
hibits three spectral features, corresponding to each ion 
addressed separately. The spectral splitting of the ions is 
given by the gradient of the coupling field, which in this 
case is proportional to the voltage applied on the RF 
electrodes. Using normal operation voltages, one reaches 
easily several kHz separation. This is an important figure 
of merit since small splittings require correspondingly 
longer pulses in order to frequency-resolve the peaks. 
Comparing to the method of individual addressing using 
magnetic field gradients @, 0|, the maximum splitting 
that we reached, of about 10 kHz, corresponds to a 
magnetic field gradient of about 50 G/cm (assuming 
a magnetic field sensitivity of 2.8 MHz/G), which is 
similar to the one obtained in [9]. The Rabi frequency 
for each ion is easily deduced from a spectrum such as 
Figj3^, taken with a different pulse time. Such a 2D scan 
is displayed in Figj3jD. We observe that the 7r-time of all 
ions are inversely proportional to their addressing fre- 
quency, since Xj — aQj, we have = ir/(aQj) (dashed 
black hyperbolas). Using the first ion frequency and its 
corresponding 7r-time as the only fitting parameters (the 
splitting for the trap parameters being deduced from an 
independent measurement as in Figj2t), we calculated 
the map showed in Figj3fc, which agrees reasonably well 
with the experimental data. For the 4 kHz splitting 
shown in Figj3^, a 7r-pulse of respectively 200, 280 and 
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FIG. 3: (Color online) Spectra of the single- addressing pro- 
tocol, (a): Spectrum for a three- ion chain taken after a 
pulse time of 400 fis, in a trap whose secular frequencies were 
(uj x ,ujy,Lj z ) = 2tt x (2.57, 2.45, 1.68) MHz. The RF frequency 
of the trap was 21.75 MHz. (b) and (c) Experimental and 
theoretical spectra taken for various driving pulse times. Ev- 
ery row in (b) is a similar spectrum to (a). The different Rabi 
frequencies for each ion are evident. The red arrows mark the 
different resonances of different ions. The color bar refers to 
the probability of having one and only one dark ion. 

560 /is leads to a spurious excitation of adjacent ions, 
which results from the overlap of their spectral response. 
The total crosstalk error, denned as the probability of 
preparing the whole chain in a different state than the 
desired one is estimated from Eq.((2]) to be respectively 
1.6,1.1 and 4.1 % when addressing each ion respectively. 
The addressing fidelities of 90(2), 88(2), 85.6(2. 5) % are 
consistent with an additional state preparation error per 
ion of 4 %, due to the spin- lock efficiency limited at long 
holding times. This is also consistent with the fidelity 
reached in a two-ion chain, of 94(2) % (with 7r-time of 
180 /is). This limitation is not fundamental, and could 
be mitigated by increasing the power in the dressing 
field, or equivalently by increasing the RF voltage, both 
of which would increase the frequency gradient, and 
allow for shorter <j x operations. 

Lastly, we verify the spatially-selective nature of this 
protocol by taking images on an Electron Multiplying 
Charge Coupled Device (EMCCD) camera. In the upper 
panel of FigHJ we show an image of a three-ion chain, 
where all the ions are in the bright state. We then 
selectively flip each ion spin using the above protocol, 
and average over 50 images. We indeed observe success- 
ful single-spin addressing. The faint "ghost" images of 
adjacent ions are due to errors (see caption of FigH]). 



FIG. 4: (Color online) Imaging the single-spin addressing pro- 
tocol on a three-ion chain. Each image below is an average of 
50 pictures taking on an EMCCD camera, with an exposure 
time of 20 ms. In the upper image, all three ions are in the 
S state, and therefore fluoresce. In the lower images a single 
ion remained in the S state while the other two were trans- 
ferred to the D (dark) state. Because of the long exposure 
time required for the EMCCD imaging, there is a significant 
amount of error due to the finite lifetime of the dark state, in 
addition to the crosstalk errors. For instance on the second 
image, the error is respectively 7(4) % and 12(5) % for the 
central and rightmost ion respectively. 



To conclude, the method described above is general, 
not limited to small systems, and relies on single-spin 
operations only. It allows the generation of any product 
state of single-qubit operations, and can in principle be 
adapted to various quantum systems. With trapped 
ions, this scheme is well suited to address single ions 
in quantum registers using the typical intrinsic in- 
homogeneity of micromotion in RF Paul traps and 
without any additional elements or any need for optical 
resolution. Our method could be useful with continuous 
density distributions as well, to select atoms with respect 
to their coupling to some external field. In particular, 
atoms with a narrow-line transition, such as neutral 
Sr or Yb, could be addressed with a sub- wavelength 
resolution, using the spatial inhomogeneity of the beam 
itself. Similarly, our technique could be readily useful in 
atomic fountain clocks to select atoms according to their 
coupling in the preparation/interrogation microwave 
cavities, where the inhomogeneities arise from the mi- 
crowave spatial mode, and would provide an additional 
knob to study collisional or cavity phase shifts in these 
clocks. 
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During the preparation of this manuscript, we became 
aware of a related work, realizing individual- ion address- 
ing of a two-qubit chain using microwave field gradients 
on a microfabricated trap [21] . 
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